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Abstract 

The QCD partition function in the external stationary gluomag- 
netic field is computed in the third order in external field invariants 
in arbitrary dimension and arbitrary covariant gauge. The contribu- 
tions proportional to third order invariants in gluon field strength are 
shown to be dependent on covariant quantum gauge fixing parameter 
a. 
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1 Introduction 



The notion of effective action plays a central role in exploring the nonlinear 
effects in quantum field theory p!j. This language especially convenient for 
analyzing the situations in which a description can be separated into contri- 
butions of quantum modes and external fields. Taking into account quantum 
fluctuations changes the tree-level description of the system through new 
nonlinear contributions in the tree-level fields that effectively correspond to 
new terms in the action of the system. Technically the computation is done 
for quantum modes propagating in the corresponding external field. These 
computations are, in particular in Yang-Mills theory, rather nontrivial. The 
physical answer should ideally be exphcitly gauge invariant with respect to 
gauge transformations of both tree-level and quantum gauge field configu- 
rations. The invariance with respect to gauge transformations of classical 
fields is ensured by using the background field gauge P^. Much more diffi- 
cult is to cope with the dependence on the quantum gauge fixing parameter. 
Generically nothing forbids a dependence of contributions proportional to 
higher order invariants in the external field being from the quantum gauge 
fixing procedure chosen. Explicit calculations of corresponding contributions 
in high energy QCD [SI El illustrate this point. An interesting suggestion of 
coping with such problems is to consider a so-called gauge invariant effective 
action j3], which in the case of Yang-Mills theory amounts to performing 
calculations in the Landau gauge a = 0. 

Of special interest is calculation of the effective action on manifolds with 
cylindrical geometry corresponding to computing the partition function (free 
energy) of the theory under consideration at finite temperature. The calcu- 
lation of the one-loop QCD partition function in external gluomagnetic field 
including the third order invariants in the external gluon field strength in 
the Feynman gauge a = 1 was performed in [H]. Fully resummed nonlocal 
contributions to a number of lowest order contributions was computed in 
[7]. Let us note that apart from having purely theoretical significance, anal- 
ysis of nonlinear contributions to the effective action (free energy) at finite 
temperature is of clear phenomenological interest. Nonlinear terms lead, in 
particular, to modification of the speed of sound in the nonabelian medium 
jH] which affects its hydrodynamical evolution [9 . 

In the present paper we generalize the calculation of |^ to the case of 
arbitrary dimension and arbitrary covariant gauge. 
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2 Computation of the partition function 



The expression for the one-loop partition function (Z) of gluodynamics in 
the covariant gauge reads 

"f^'^ZM = -J^/<i'"2^(Gy'-5lndetW(a) + hidet(--^C=^) 

- (1) 

Here F is free energy, T is the temperature, is a kinetic operator for the 
ghosts, and is a gluon kinetic operator 

= -D\Ar5,. - 2r''^G% - - l)D;^Df (2) 

where a is a gauge fixing parameter, /'^'"^ are the structure constants of 
the corresponding Lie algebra, is a covariant derivative containing the 
external gauge field potential and G^j, is the corresponding field strength. 
Let us emphasize that we keep in ((H) a manifest dependence on the gauge 
parameter a-dependence in the ghost determinant (usually it is hidden in the 
normalization constant) because of its importance in performing calculations 
in the Landau gauge a — > 0. 

In computing the logarithms of the determinants in it is convenient 
to use a C - regularization 



hi det A = - lim — — - / dt t^-^Tre"*^ (3) 

As the closed expression for the right-hand side of (jHl) does not exist, one is 
forced to use an expansion of Tre~*^ in the invariants of the background field 
(Seeley expansion) 10,. For example, at zero temperature the expansion of 
gluon kinetic operator W reads 



/oo 
d^^x) w, a\x)s-^^^ (4) 

fc=0 

To calculate the functional trace in it is convenient to use the basis of 
plane waves (see, e.g., [T2].[T! 



Tr < x\e-^'\x > = Tr / — e'^^'" < x\e-^'\x > e'^" 

(27r)2'^ I I 



Tr J d'^^x J ^|^e-^^(^---^-+*P-)l , (5) 
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where the last trace is performed over the Lorenz and color indices. The ex- 
plicit expressions for the first four Seeley coefficients at zero temperature were 
computed in [HI Unj ( we refer for the detailed description of the calculation 
method to these papers) 

b-^ = f:^(2--[l-«1), 

b-^^^ = ^{2-^(2--[l-«--^])}G,. (6) 



b-uj+3 

where 



^{2uj - [1 - a--'']){Gs - 3/3) + + ^oo, a)h 



/^a r^a pabc/^a /^b r^c j ( -r^ab/^b \( -r^ac/^c \ ( '-7\ 

are the invariants one can construct from gluon field strength tensor and 

e(uj,a) = — ^ -I - \Auj-l-a ^[1 - a'^^^l ) 

' ^ 4(tu- 1) \cu V 1 -tt 7 

3a 1 - a^-^ 1 1 - a^-"^ 



l-a uj-2 J 6(cj-2) 

represents nontrivial off-shell contribution vanishing in the Feynman gauge 
at a = 1. 

To calculate the trace in ^ at the finite temperature it is convenient to 
use the Matzubara formalism in which the Eucledian space has cylindrical 
geometry. The extension of the compact dimension is the inverse temperature 
(3 = 1/T and the fields are periodic in this coordinate with the period j3. As 
a result the temporal integral in is replaced by the sum over Matsubara 
frequencies ujn = {211/ [3)n and © takes the following form 

Since the case of zero Matsubara frequency n = is completely analogous to 
the zero temperature case for the theory in dimension 2^1; — 1 in what follows 
we consider only the temperature-dependent contributions I]^_oo ^ 2 Z^J^o 
and keep only even powers of Un- 
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Technically summation in (jH)) can be performed in closed form only to- 
gether with integration over proper time in Q. 

Let us first consider the contributions to the free energy in the first non- 
trivial order D'^. For the gluonic contribution one has^ 



1, d T(s + 1/2 + 2 -uj)C(2s + 5- 2uj) 
hm 



X 



2 s^o ds 
[ 47r2T2 



r s) 



(9) 



2 - 



2iu-l + a' 
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The ghost contribution in the same order Z)^ reads 

d V{s + 1/2 + 2- uj) C(2s + h-2uj) 



-d: 



lim 



X 



-0 ds 



V{s) 



(10) 



The contributions to the free energy of gluons and ghosts in the next order 
read (here when taking a limit s — we assume that 7^ 3) 



f): 



C(7-2cj)r(l/2 + 3-c.^) 



X 



T7-2^ 1180 



^'':(G3-3/3) + ^/3 + (i 



:ii^ 



7 + « 
^ 24 ^ 



and 



-F\ _ C(7 - 2cj)r(l/2 + 3 - cu) ^=^iV, 
~ ' f)& ~ 2H^-^J1^ T7-2^ il80 



3/3 



(12) 



Combining the gluon and ghost contributions for the free energy Eqs. (jn|l- (fT^ 
we obtain the final expression for the regularized free energy 



TJ 



p\9h 



TJ 



d T{s + 1/2 + 2 -uj)C{2s + 5 -2uj) 
hm- — — : — ;r^^ (13) 



^0 ds 



1677' 



A-uj 



X 



2iu 



Tis) 
a' - 2a^/2 



12 



V47r2T2 J 
C(7-2cu)r(l/2 + 3-cj) g^Nc 



Go 



2%6 

'2(lj - 1 



X 



180 



(G'3-3/3) + ^/3 + (l-«)^/3 



""^In what follows we assume that invariants ^ include integration over corresponding 
spatial coordinates 
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We see that in the last formula the a and u dependences are factorized 
unlike the zero temperature case 

Special care has to be exercised when computing the QCD partition func- 
tion in the physical dimension d = 4 {u = 2). Let us first consider the con- 
tributions of the fourth order in covariant derivatives. The expressions for 
the gluon and ghost contribution to the free energy read 

- - =^^— lOn— --na — G2 14 
VTA 167r2 3 V T2 2 J AT ^ ^ ' 

and 

- - =^^— n — + nVa — G2 . 15 
Vr/4 167r2 3 V T2 ^ y 4T ^ ' 

We see that the expressions Eqs. fll4j) and ()15|1 contain contributions propor- 
tional to the logarithm of the quantum gauge fixing parameter logo; which 
are singular in the limit a — 0. However, these terms in Eqs. ()14j) and (fT3j) 
cancel each other upon adding the gluon and ghost contributions. There are 
no other contributions singular in a in the next order in covariant derivatives 
in any dimentions^. Let us emphasize again that all singularities in a disap- 
pear only after regularization and taking into account the ghost contribution. 

At zero temperature we have the leading singularites of the form [1 — 
q;'^"^] and "subleading" terms which come from the contribution of order 

in Eq. (jH} that in the limit of u —* 2 take the following form: 

e(2, a) = ^{3(7 -a) + 2^^^(2 + 7a)}. (16) 

This basic difference between the zero and finite temperature cases is due 
to the existence of the natural energy scale (T) in the last case. If we try 
to estimate the corresponding determinants following ^3] by introducing the 
energy scale = —60/^1 and writing 

IndetW = -6o(ln^ + 7i<; - 1) (17) 

{'Je is Euler's constant) we will have the same cancellation of a ^ singular- 
ity. To elucidate this statement let us consider a variation of the functional 
determinant under scale transformation. By using ^-function definition of a 
determinant it is easy to get the following equation P: 

^At a; = 3 we have no singularities in a in the order D* and the same logarithmic 
cancellation in the order as for w = 2 in the order Z?* 
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lndet(-5) =lndet(5) -6olna. (18) 
a 

If now we estimate both of the above determinants by using the formula 
© directly we obtain, instead of (fTS|) . the following expression: 

oo 

5: (1 - a--^)A'^-^6_^+;v ~ -ho In a 

N=l 

where 1/A is the (upper) cutoff parameter in the integral over proper time. 
By comparing this expression with the explicit form of the Seeley coefficients 
© we conclude that all [1 — a"~^]-terms give nothing else than the contri- 
bution —Bq In a in the final answer. 

Summing all contributions to the free energy we obtain the final expres- 
sion for the logarithm of the QCD partition function in the physical dimen- 
sion d = 4 (cu = 2): 

, Z\A,T] ( , 1 

n \ ' \ = - 1 - ^^ — 11 n — — G2 19 
Z[0,T] \ 167r2 3 AT ^ ^ ' 



Let us stress that although the formula ()19|) does not contain singularities 
in the quantum gauge fixing parameter a, the answer does depend on it off- 
shell. 

Expressions in Eqs. (fT^ and (fTIH) constitute the main results of the paper. 



3 Conclusion and outlook 

Let us formulate the main result of the present paper. Using the heat kernel 
expansion we have computed the partition function of pure Yang- Mills theory 
to the third order in the invariants in the external field strength in arbitrary 
dimension and arbitrary covariant gauge in the case when only stationary 
spatial ("magnetic") gauge potentials of the external field were taken into 
account. The calculation demonstrates delicate a - dependent cancellations 
between the gluon and ghost contributions to the partition function. The 
first nontrivial nonlinear contribution to the effective action was shown to 
depend on the quantum gauge fixing parameter a. 

The major question left unanswered at this stage is to complete the pro- 
gram by calculating the partition function taking into account the temporal 
gauge connection Aq. Work in this direction is currently in progress [16j . 
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